CERN-TH/2001-018 
KA-TP-2-2001 
February 1, 2008 
hep-ph/0101260 



Two-loop electroweak contributions to Ar Q 



A. Freitas", S. Heinemeyer 6 , W. Hollik c , W. Walter c , G. Weiglein^ 

a DESY Theorie, Notkestr. 85, D-22603 Hamburg, Germany 

b HET, Physics Department, Brookhaven Nat. Lab., Upton, NY 11973, USA 

c Institut fur Theoretische Physik, Universitat Karlsruhe, 
D-76128 Karlsruhe, Germany 

d CERN, TH Division, CH-1211 Geneva 23, Switzerland 



A review is given on the quantum correction Ar in the W-Z mass 
correlation at the electroweak two-loop level, as derived from the calcula- 
tion of the muon lifetime in the Standard Model. Exact results for Ar and 
the PF-mass prediction including 0(a 2 ) corrections with fermion loops are 
presented and compared with previous results of a next-to-leading order 
expansion in the top-quark mass. 



Presented by W. Hollik at the 

5th International Symposium on Radiative Corrections 

(RADCOR-2000) 
Carmel CA, USA, 11-15 September, 2000 



*Work supported by the Deutsche Forschungsgemeinschaft and by the European Union. 



1 Introduction 



The interdependence between the W^-boson mass, Mw, and the Z-boson mass, 
Mz, with the help of the Fermi constant Gf and the fine structure constant a is 
one of the most important relations for testing the electroweak Standard Model (SM) 
with high precision. At present, the world-average for the M^-boson mass is M^ p = 
80.434 ± 0.037 GeV [jlj. The experimental precision on My/ will be further improved 
with the data taken at LEP2 in their final analysis, at the upgraded Tevatron and at 
the LHC, where an error of 5Mw = 15 MeV can be expected @. At a high-luminosity 
linear collider running in a low-energy mode at the W + W~ threshold, a reduction 
of the experimental error down to 5My/ = 6 MeV may be feasible 0. This offers 
the prospect for highly sensitive tests of the electroweak theory jl| , provided that the 
accuracy of the theoretical prediction matches the experimental precision. The basic 
physical quantity for the My/-Mz correlation is the muon lifetime r M , which defines 
the Fermi constant Gf according to 



1 G 



2 ,^5 




1 + A QED ), (1) 



with F(x) = 1 — 8x — 12x 2 In x + 8x 3 — x A . By convention, the QED corrections within 
the Fermi Model, Aq ED , are included in this defining equation for Gp- The one-loop 
result for Aqed [§], which has already been known for several decades, has recently 
been supplemented by the two-loop correction |J, yielding 

A QED = 1 - 1.81 « + 6.7 f ^>V . with a(m „) * _L- . (2) 

The tree- level jy-propagator effect giving rise to the (numerically insignificant) term 
3m^/(5M^/) in (P, is conventionally also included in the definition of Gp, although 
not part of the Fermi Model prediction. From the precisely measured muon-decay 
width the value | G F = (1.16637 ± 0.00001) 10~ 5 GeV" 2 for the Fermi constant is 
derived. 

Calculating the muon lifetime within the SM and comparing the SM result with (0) 
yields the relation 

where the radiative corrections are summarized in the quantity Ar, first calculated 
in H at one-loop. This relation can be used for deriving the prediction of Mw within 
the SM (or possible extensions), to be confronted with the experimental result for 
M w . 
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The one-loop result for Ar within the SM can be decomposed as follows (with the 
notation 4 = 1- M^/Mf , = 1 - ^ 



_w 

-.2 
>w 



Ar^ = Aa - ^Ap + Ar rcm (M H ), (4) 



exhibiting the leading fermion-loop contributions Aa and Ap, which originate from 
the charge and mixing-angle renormalization; the remainder part Ar rem contains in 
particular the dependence on the Higgs-boson mass, Mjj. The QED-induced shift Aa 
in the fine structure constant contains large logarithms of light-fermion masses. The 
leading contribution to the p parameter from the top/bottom weak-isospin doublet, 
Ap, gives rise to a term with a quadratic dependence on the top-quark mass, m t |J . 

Beyond the one-loop order, resummations of the leading one-loop contributions 
Aa and Ap are known They correctly take into account the terms of the form 
(Aa) 2 , (Ap) 2 , (AaAp), and (AaAr rem ) at the two-loop level and the leading powers 
in Aa to all orders. 

QCD corrections to Ar are known at 0(aa s ) [|TTJ and 0(aa 2 s ) |12|. Concerning the 



electroweak two-loop contributions, only partial results are available up to now. Ap- 
proximative calculations were performed based on expansions for asymptotically large 
values of M H JTB) and m t [Tjrr^|T^| . The terms derived by expanding in the top-quark 



mass of 0(G F rn,t) [0 and 0(G F mfM z ) JT5| were found to be numerically sizeable. 
The 0(G F m 2 M^) term, involving three different mass scales, has been obtained by 
two separate expansions in the regions M w , M z , M H -C m t and M w , M z <C m t , M H 
and by an interpolation between the two expansions. This formally next-to-leading 
order term turned out to be of a magnitude similar to that of the formally leading 
term of 0(G 2 F mf), entering with the same sign. Its inclusion (both for M\y and the 
effective weak mixing angle) had important consequences on the indirect constraints 
on the Higgs-boson mass derived from the SM fit to the precision data. 

A more complete calculation of electroweak two-loop effects is hence desirable. As 
a first step in this direction, exact results were derived for the Higgs-mass dependence 
of the fermionic two- loop corrections to the precision observables Wh - They have been 
compared with the results of expanding up to 0{G F m 2 M z ) |T5|| , specifically analysing 
the effects of the m t expansion, and good agreement has been found [|18|]. Beyond 
the two-loop order, complete results for the pure fermion-loop corrections {i.e. contri- 
butions containing n fermion loops at ra-loop order) have recently been obtained up 
to four-loop order |I9 |. These results contain in particular the contributions of the 



leading powers in Aa as well as the ones in Ap and the mixed terms. 

In this talk the complete fermionic electroweak two-loop corrections to Ar are 
discussed, as calculated exactly without an expansion in the top-quark or the Higgs- 
boson mass P0|. These are all two-loop diagrams contributing to the muon-decay 
amplitude and containing at least one closed fermion loop (except the pure QED 
corrections already contained in the Fermi model result, according to ([[])). Fig. |l| 
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displays some typical examples. The considered class of diagrams includes the po- 
tentially large corrections both from the top/bottom doublet and from contributions 
proportional to Nif and Nff, where Nif is the number of light fermions (a partial 
result for the light-fermion contributions is given in [21]). The results presented here 
improve on the previous results of an expansion in mt up to next-to-leading order |j| 
in containing the full dependence on m t as well as the complete light-fermion contribu- 
tions at the two-loop order, while in [HJ higher-order corrections from light fermions 
have only been taken into account via a resummation of the one-loop light-fermion 
contribution. 




2 Outline of the calculation 

Since all possibly infrared (IR) divergent photonic corrections are already con- 
tained in the definition of the Fermi constant Gf and mass singularities are 
absorbed in the running of the electromagnetic coupling, My/ represents the scale 
for the electroweak corrections in Ar. Therefore it is possible to neglect all fermion 
masses except the top-quark mass and the momenta of the external leptons so that 
the Feynman diagrams for muon decay reduce to vacuum diagrams. 

All QED contributions to the Fermi Model have to be excluded in the computation 
of Ar since they have already been separated off in the definition of Gp, see eq. ([!]). 
Apart from the one-loop contributions, this comprises two-loop QED corrections and 
mixed contributions of QED and weak corrections of each one-loop order, which 
have to be removed from Ar^ Q For fermionic two- loop diagrams it is possible 
to find a one-to-one correspondence between QED graphs in Fermi-Model and SM 
contributions. 

After extracting the IR-divergent QED corrections, the generic diagrams con- 
tributing to the muon-decay amplitude can be reduced to vacuum-type diagrams, 
since the masses of the external particles and the momentum transfer are negligi- 
ble. The renormalization of is performed in the on-shell scheme. Thus, the mass- 
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renormalization of the gauge bosons requires the evaluation of two-loop two-point 
functions with non-zero external momentum, which is more involved from a technical 
point of view regarding the tensor structure and the evaluation of the scalar integrals. 
This complication cannot be avoided by performing the calculation within another 
renormalization scheme, e.g. the MS scheme, since ultimately one is interested in 
the relation between the physical parameters Mw, Mz, a, Gp, where the two-point 
functions for non-zero momenta enter. 

All diagrams and amplitudes for the decay and counterterm contributions have 
been generated with the program FeynArts 2.2 |[22|| . The amplitudes are algebraically 
reduced by means of a general tensor-integral decomposition for two-loop two-point 
functions with the program TwoCalc [53], leading to a fixed set of standard scalar 



integrals. Analytic expressions are known for the scalar one-loop |24| and two-loop 
75| vacuum integrals, whereas the two-loop self-energy diagrams can be evaluated 



numerically by means of one- dimensional integral representations [26 



In order to apply an additional check the calculations were performed within a 
covariant gauge, with individual gauge parameters & for each gauge boson. It has 
been explicitly checked at the algebraic level that the gauge-parameter dependence 
of the final result drops out. 

At the subloop level, also the Faddeev-Popov ghost sector has to be renormalized. 
The gauge- fixing part of the Lagrangian, in terms of the gauge fields A^, Z^, W 
and the unphysical Higgs scalars Xi ^ given by 



C 



gf 



(F 7 ) 2 + (F z f + F + F~ + F~F^ 



with 



F^={CiY h d^ + ^d,Z^ 



(5) 



F Z = (ef ) 1 d»Z» + ^—d^ - (H) 1 M z X , 

F ± = (er)"' d^W* =F i h M w ^, 

does not need renormalization. Accordingly, one can either introduce the gauge-fixing 
term after renormalization or renormalize the gauge parameters in such a way that 
they compensate the renormalization of the fields and masses. Both ensure that no 
counterterms arise from the gauge-fixing sector but they differ in the treatment of 
the ghost Lagrangian, which is given by the variation of the functionals F a under 
infinitesimal gauge transformations 58b, 

X pa 

o,o=7,Z,± 



4 



In the latter case, which was applied in our work for simplification of the automatized 
treatment, additional counterterm contributions for the ghost sector arise from the 
gauge-parameter renormalization. The parameters in @ are renormalized such 
that their counterterms 5^1 exactly cancel the contributions from the renormalization 
of the fields and masses and that the renormalized gauge parameters comply with the 
gauge. 



3 On the 75— problem 

In four dimensions the algebra of the 75-matrix is defined by the two relations 

{75,7j = for a = l,...,4 (7) 

Tr{757'V7V} = Aie^ pa . (8) 

It is impossible to translate both relations simultaneously into D 7^ 4 dimensions 
without encountering inconsistencies f27j. 



A certain treatment of 75 might break symmetries, i.e. violate Slavnov- Taylor 
(ST) identities which would have to be restored with extra counterterms. Even after 
this procedure a residual scheme dependence can persist which is associated with 
e-tensor expressions originating from the treatment of (H). Such expressions cannot 
be canceled by counterterms. If they broke ST identities this would give rise to 
anomalies. 



't Hooft and Veltman [27| suggested a consistent scheme, formalized by Breiten- 
lohner and Maison |28j , as a separation of the first four and the remaining dimensions 
of the 7-Matrices (HVBM-scheme). It has been shown that the SM with HVBM 
regularization is anomaly-free and renormalizable. This shows that e-tensor terms do 
not get merged with divergences. 

The naively anti-commuting scheme, which is widely used for one-loop calcula- 
tions, extends the rule (^) to D dimensions but abandons 

{75,7a} = for a = l,...,D (9) 

Tl { l5l »YYr} = 0. (10) 

This scheme is unambiguous but does not reproduce the four-dimensional case. 

In the SM, particularly triangle diagrams (like the ones in Figure ^) containing 
chiral couplings are sensitive to the 75-problem. In our context, the one- loop triangle 
diagrams have been explicitly calculated in both schemes. While the naive scheme 
immediately respects all ST identities the HVBM scheme requires the introduction 
of additional finite counterterms. Even after this procedure finite differences remain 
between the results of the two schemes, showing that the naive scheme is inapplicable 
in this case. 
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Figure 2: Charged-current vertex diagrams with fermion-triangle subgraphs. 



In the calculation of Ar triangle diagrams appear as subloops of two-loop charged 
current (CC) vertex diagrams (Fig. One finds that for the difference terms be- 
tween both schemes this loop can be evaluated in four dimensions without further 
difficulties. This can be explained by the fact that renormalizability forbids diver- 
gent contributions to e-tensor terms from higher loops in the HVBM scheme. The 
e-tensor contributions from the triangle subgraph in the HVBM scheme meet a second 
e-tensor term from the outer fermion lines in Fig. [| thereby resulting in a non-zero 
contribution to Ar. 

Computations in the HVBM scheme can in general get very tedious because of the 
necessity of additional counterterms. For our specific problem, however, it is possible 
to apply another, simplified, method. One can consider a "mixed" scheme that uses 
both relations (0) and (]8|) in D dimensions, despite their mathematical inconsistency, 
to evaluate the one-loop triangle subgraphs. These results immediately respect all 
ST identities. As checked explicitly, they differ from the HVBM results (with the 
appropriate counterterms to restore the ST identities) only by terms of 0(D — 4), 

r HVBM = rm ix ) + 0(jD _ 4) _ (n) 

Inserting the one-loop expressions into the the two-loop diagrams one finds that the 
second loop integration gives a finite result and hence can be performed in four 
dimensions yielding 

r^ = r cc (2) + 0p-4). (12) 

Thus the mixed scheme can serve in this case as a technically easy prescription for 
the correct calculation of the CC two-loop contributions. Practical ways of treating 
75 in higher-order calculations are also discussed in |5D| . 



4 Two-loop renormalization 



For the determination of the one-loop counterterms and renormalization constants 
the conventions of |3l| are adopted. Two-loop renormalization constants enter via the 
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counterterms for the transverse W propagator and the charged current vertex (the 
counterterms for the transverse Z propagator are analogous): 



w w 



5Z*(k 2 



M 2 W ) 



(13) 




V2& 



1^- 



sz, 



e(2) 



^W(2) 



i/h 



+ (1-loop renormalization constants) 



(14) 



SZ w,eL,uL denote the field-renormalization constants, 5M^ Z the W- and Z-mass coun- 
terterms, and SZ e denotes the charge-renormalization constant. The lower indices in 
parentheses indicate the loop order. The mixing-angle counterterm 5sw(2) can be 
derived from the gauge-boson mass counterterms. The two-loop contributions always 
include the subloop renormalization. 

The on-shell masses are defined as the position of the propagator poles. Starting at 
the two-loop level, it has to be taken into account that there is a difference between the 
definition of the mass M 2 as the pole of the real part of the (transverse) propagator, 



Re{( J D T )- 1 (Af 2 )} = 0, 

and the real part M 2 of the complex pole, 

{D T y\M 2 ) = 0, M 2 = M 2 -iMT. 



(15) 



(16) 



The imaginary part of the complex pole is associated with the width V. The defining 
condition ( |i~6l) yields for the VT-mass counterterm 



5M, 



W(2) 



Re{^ {2) (M 2 w )}-5Z^SM^ {1) 



+Im{E^ / 1) (M^)}Im{S^ 1) (M^)}, (17) 



whereas for the real-pole definition the last term of eq. (|17] ) is missing. denotes 
the transverse W self-energy and its momentum derivative. Similar expressions 
hold for the Z boson. 

The W and Z mass countertermss determine the two-loop counterterm for the 
mixing angle, <5sw(2), which has to be gauge invariant since s w is an observable quan- 
tity. With the use of a general gauge it has explicitly been checked that 5sw(2) 
is gauge-parameter independent for the complex-pole mass definition, whereas the 
real-pole definition leads to a gauge dependent <5sw(2)- This is in accordance with the 
expectation from S'-matrix theory [fj^j , where the complex pole represents a gauge- 
invariant mass definition. 

It should be noted that the mass definition via the complex pole corresponds to 
a Breit-Wigner parameterization of the resonance shape with a constant width. For 
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the experimental determination of the gauge-boson masses, however, a Breit-Wigner 
ansatz with a running width is used. This has to be accounted for by a shift of the 
values for the complex pole masses 



_ r 2 

M = M - 



2M y ' 

which yields the relations 

M_ z = M z - 34.1 MeV, . , 

M w = M w - 27 A (27.0) MeV for M w = 80.4 (80.2) GeV. 1 } 

For Mz and V z the experimental numbers are taken. The W mass is a calculated 
quantity, and therefore also a theoretical value for the H^-boson width should be 
applied here. The results above are obtained from the approximate, but sufficiently 
accurate expression for the W width, 



w 



,G F M W , 2a, 



2V2tc V 3tt 



1 + . (20) 



5 Results 

In the previous sections the characteristics of the calculation of electroweak two- 
loop contributions to Ar have been pointed out. Combining the fermionic 0(a 2 ) 
contributions with the one-loop and the QCD corrections yields the total result 

Ar = Ar (a) + Ar K) + Ar {aa ^ + Ar (N f a2) + Ar (JV /° 2) . (21) 

Here Nf,Nj symbolize one and two fermionic loops, respectively. Fig. |^ shows that 
both the QCD and electroweak two-loop corrections give sizeable contributions of 
10-15% with respect to the one- loop result. 

In Fig. § the prediction for Mw derived from the result (|2l|) and the relation 
(0) is compared with the experimental value for M w . Dotted lines indicate one 
standard deviation bounds. The main uncertainties of the prediction originate from 
the experimental errors of m t = (174.3 ± 5.1) GeV 0] and Aa = 0.05954 ± 0.00065 
f34|. It is obvious that low Higgs masses are favored; the new results on Ar strengthen 
the tendency towards a lighter Higgs boson (according to the following comparison). 

These results can be compared with the results obtained by expansion of the two- 
loop contributions up to next-to-leading order in m t flT5j. The predicted values for 



M\y for several values of M# are given in Tab. [I]. For the input parameters the values 
of @ have been chosen, i.e. m t = 175 GeV, M z = 91.1863 GeV, Aa = 0.0594, 
a s (Mz) = 0.118. Agreement is found between the results with maximal deviations of 
less than 5 MeV in M^y. The deviations in the last column of Tab. [I] can of course not 
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Figure 3: Various stages of Ar, as a function of Mh- The one-loop contribution, AA a \ is 
supplemented by the two-loop and three-loop QCD corrections, Arq^, D = Ar' ms )+Ar' aa ^, 
and the fermionic electroweak two-loop contributions, Ar^" 2 ) = AA N{a2 ^ + Ar^f" 2 '. For 
comparison, the effect of the two-loop corrections induced by a resummation of Aa, Ar^" a , 
is shown separately. 



be attributed exclusively to differences in the two-loop top-quark and light-fermion 
contributions, because the results also differ by a slightly different treatment of those 
higher-order terms that are not yet under control, such as purely bosonic two- loop 
contributions, and effects from scheme dependence. 

Similar to |35||, a simple numerical parametrization of our result for Mw can be 



given by the following expression: 

M w = M w - ci dH - c 5 dH 2 + c 6 dH 4 -c 2 da + c 3 dt - c 7 dH dt - c 4 da s , (22) 

where 

dH=m(^-v dt= r m < 

V 100 GeV/' V 174.3 GeV 

da = -^ 7 -l, (23) 
0.05924 ' 0.119 ' V 7 

with the coefficients M<b, = 80.3767 GeV, c x = 0.05613, c 2 = 1.081, c 3 = 0.5235, 
c 4 = 0.0763, c 5 = 0.00936, c 6 = 0.000546, c 7 = 0.00573. and with M z = 91.1785 



GeV. The quality of the approximation (22) to our full result for Mw is within 



0.4 MeV, allowing M H between 65 GeV and 1 TeV. 
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80.5 





M™ p = (80.434 ± 0.037) GeV ■ 
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M H [GeV] 

Figure 4: The SM prediction for Mw as a function of Mh for mt = 174.3 ±5.1 GeV is 
compared with the current experimental value, M^ p = 80.434 ± 0.037 GeV [jl]. 

6 Conclusion 

In this talk the realization of an exact two-loop calculation of fermionic contribu- 
tions in the full electroweak SM and its application to the precise computation of Ar 
has been reviewed. Numerical illustrations were given for the results, which might 
serve as ingredients for future SM fits. 



Table 1: Comparison between M^y-predictions from an NLO expansion in mt (M^ pa ) and 
the full calculation (M^ n ). 5Mw denotes the difference. 



M H 


M^ pa 


M^ 11 


5M W 


[GeV] 


[GeV] 


[GeV] 


[MeV] 


65 


80.4039 


80.3997 


4.2 


100 


80.3805 


80.3771 


3.4 


300 


80.3061 


80.3051 


1.0 


600 


80.2521 


80.2521 


0.0 


1000 


80.2129 


80.2134 


-0.5 
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